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Abstract. - Periodic orbit quantization requires an analytic continuation of non-convergent 
semiclassical trace formulae. We propose a method for semiclassical quantization based upon the 
Pade approximant to the periodic orbit sums. The Pade approximant allows the re-summation 
of the typically exponentially divergent periodic orbit terms. The technique does not depend 
on the existence of a symbolic dynamics and can be applied to both bound and open systems. 
Numerical results are presented for two different systems with chaotic and regular classical 
dynamics, viz. the three-disk scattering system and the circle billiard. 
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Semiclassical theories link the spectrum of a quantum system to the dynamics of its classical 
counterpart and thereby play an important role in the deeper understanding of the relation 
between quantum and classical mechanics. Of particular interest are periodic orbit theories 
which describe the quantum mechanical density of states in terms of contributions from the 
periodic orbits of the classical system. The semiclassical trace formulae have been derived by 
Gutzwiller for classically chaotic systems [[j} |[ and by Berry and Tabor for regular (integrable) 
systems (j| . A common feature of these formulae is that in most cases the periodic orbit sum 
does not converge in those regions where the semiclassical eigenenergies or resonances are 
located. These resonances are given as the poles of the periodic orbit sum, 

g(E) =Y / A po(E)e l[s o° {E)/R -%^° ] , (1) 

po 

where A po , S po , and fi po are the amplitude, classical action and Maslov index of the periodic 
orbit (po), respectively. Formulae for the amplitudes A po are given in refs. |^|, H|. 

In the last decade various techniques have been developed to circumvent the convergence 
problem of periodic orbit theory. Examples are the cycle expansion technique ^, ||, |^, |7) , the 
Riemann-Siegel type formula and pseudo-orbit expansions pi , and a periodic orbit quantization 
rule for the counting function of eigenvalues ^ . These techniques have been proven to be very 
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efficient for systems with special properties, e.g., the cycle expansion for hyperbolic systems 
with an existing symbolic dynamics, and the other methods for the calculation of bound 



spectra. Recently, a method for periodic orbit quantization has been introduced ]10| , 11 
[l2| which is based on the high resolution analysis (harmonic inversion) of the semiclassical 
recurrence signal. The method requires the knowledge of the periodic orbits up to a given 
maximum period (classical action), which depends on the mean density of states (i.e., the 
Heisenberg time), and can be applied to open or bound systems with regular, chaotic, or even 
mixed fl3| classical dynamics. 

In this Letter we propose an alternative method for periodic orbit quantization. It is 
the Pade approximant (PA) to slowly convergent and/or divergent periodic orbit sums. In 
the former or the latter case, the PA either significantly increases the convergence rate or 
analytically continues the exponentially divergent series. The PA is especially robust for 
re-summing diverging series in many applications in mathematics and theoretical physics 
| jl4| . An important example is the summation of the divergent Rayleigh-Schrodinger quantum 
mechanical perturbation series for, e.g., atoms in electric |l5| and magnetic fl6[| fields. In 
periodic orbit theory the PA has been applied to cycle-expanded Euler products and dynamical 
zeta functions ||. However, to the best of our knowledge, the PA has not yet been used for 
the direct summation of periodic orbit terms. 

Pade approximant (PA) to the periodic orbit sums. - The PA to a complex function f(z) 
is defined as a ratio of two polynomials and can be computed from the coefficients a n of 
the Maclaurin expansion of f(z), i.e., a power series f(z) = X^^Lo " 2 ™ W1 th finite or even 
zero radius of convergence in z |16| . However, eq. (Q) does not have the functional form of a 
Maclaurin power series expansion of g(E) in the energy E. Even disregarding this limitation, 
a direct computation of the PA to the sum in eq. (0) would be numerically unstable due to the 
typically large number of periodic orbit terms. Nevertheless, considering J? as a parameter, 
g{E) can be rearranged and written as a formal power series in an auxiliary variable z, 



£ A po (E)e lS ^/ h 



= £o»(£)* n . (2) 



where the maximal value of n required for convergence of the PA is relatively small compared 
with the number of periodic orbit terms. Of course, the arrangement (||) of the periodic orbit 
sum as a power series is not unique, and expansions similar to eq. (0) can be used for other 
ordering parameters n of the orbits, e.g., the cycle length in systems with an existing symbolic 
code. However, if no symbolic dynamics exists, the sorting of orbits by their Maslov index 
is natural both physically and as a way to introduce an integer summation index and will 
be justified below by the successful numerical application of the method to a regular system 
without symbolic dynamics. (For a discussion of the ordering of orbits in cycle-expansion 
techniques see Mj.) Note that the PA to the periodic orbit sum cannot be applied without 
any ordering parameter, i.e., when no symbolic dynamics exists and all Maslov indices are zero, 
which is the case, e.g., for the Riemann zeta function as a mathematical model for periodic 
orbit quantization p| . The true value g(E) of the periodic orbit sum is obtained by setting 
z = 1 in eq. (||), i.e., g(E) — g(l;E). In such a case, we have a point PA which is given as 
a ratio of two polynomials in z whose coefficients are non-polynomial functions of E all at 
a fixed value of z. The usual implementation of the PA as a ratio of two polynomials in z 
whose coefficients are computed, e.g., via the Longman algorithm ]l8| would be advantageous 
if g[z; E) were required for many values of z. In the present case, at each given energy E 
only one fixed value z = 1 is needed and the PA is most efficiently computed by means of the 
recursive Wynn e-algorithm jl9|. 
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To briefly describe the e-algorithm, we introduce a sequence of partial sums {A„} which 
converges to (or diverges from) its limit A as n — > oo. In the case of divergence, A is called 
the 'anti-limit' of {A n } , as n — ► oo. Further, let F be a transformation which maps {-An} 
into another sequence {B n }. The mapping F will represent an accelerator, i.e., sequence 
{£>„} will converge to the same limit A faster than {A n } if the following condition is fulfilled: 
(B n — B)/(A n — A) — > 0, as n — > oo. In addition, the same F can be applied to wildly 
divergent sequences {j4 n }. Only non-linear mappings can simultaneously accomplish both 
goals to accelerate slowly convergent and induce convergence into divergent sequences. The 
transformation F will be non-linear if its coefficients depend on A n , e.g., the so-called e- 
algorithm of Shanks (20), whose mapping F is the operator eu which converts sequence {A n } 
into {B n } via ek{A n ) = B^^ n = [to + k/k] (n > 0, n > k). This is the well-known Aitken 
A 2 -iteration process {i.e., the simplest PA, [1/1]) extended to higher orders k. The general 
term in the fcth-order transform of A n can be computed efficiently from the stable and 

recursive e-algorithm of Wynn i.e., e s (j4 m ) = = [m + s/s], where 

effi = ei m + 1} + V - 4 m) ) ; rn, s > (3) 

with — 0, £g"^ = A m , £2s+i = ^/ e s(AA m ) and where A is the forward difference 

operator: Axj = Xj +i — Xj. When {A n } is the sequence of partial sums of a power series, 

the two-dimensional array £^~ s ^ yields the upper half of the well known Pade table [m/s]. 
However, the e-algorithm need not necessarily be limited to power series. 

The procedure to apply the above PA to semiclassical quantization by summation of periodic 
orbit terms is as follows. For a given system we calculate the periodic orbits up to a chosen 
maximum ordering parameter n < n max where n can be but is not necessarily related to the 
Maslov indices of orbits. Note that this set of orbits usually differs from the set of orbits with 
classical action Sp Q < S max , which is required for periodic orbit quantization by harmonic 
inversion [jl0| |ll|, From the quantities A po , S po and /x po of these orbits, we compute the 
partial sums (e.g., with n < ^i m ax the Maslov index): 

An = A po (E)e t[s p° {E)/n -%^° ] . (4) 

The sequence {A n } of partial sums is used as input to the e-algorithm, eq. (|J), to obtain 
a converged value g(E) of the periodic orbit sum (|l|). The semiclassical eigenenergies or 
resonances are given as the poles of g(E) and are obtained by searching numerically for the 
zeroes of the reciprocal function l/g(E). Such a search requires the evaluation of, e.g., A po (E) 
and S po (E) at complex values of E. This is straightforward for the scaling systems considered 
in this Letter. The number of the poles of g(E) is not constrained by the size of the sequence 
{A n } of partial sums since our PA is a ratio of two non-polynomial functions of E. We now 
demonstrate the power of the method on two physical examples with completely different 
dynamical properties, viz. the open three-disk billiard and the bound circle billiard. 

The three-disk scattering system. - As the first example we consider a billiard system 
consisting of three identical hard disks with unit radii, R = 1, displaced from each other 
by the same distance d. This simple, albeit non-trivial, scattering system has served as 
a model for the development of the cycle expansion method [^, ||, ^, |7j] and periodic orbit 
quantization by harmonic inversion [ p)[ pT[ O]. The three-disk scattering system is invariant 
under the symmetry operations of the group Cs v , i.e., three reflections at symmetry lines and 
two rotations by 27r/3 and 47r/3. Resonances belong to one of the three irreducible subspaces 
Ai, A 2 , and E [pTI . In the following we concentrate on the resonances of the subspace A\. For 
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Fig. 1. - Convergence vs. exponential divergence of the partial periodic orbit sums for the three-disk 
scattering system with R = 1, d = 6 as functions of the order n at (a) complex wave number 
k — 150 — O.li and (b) k — 150 — 0.5i. Dashed lines and plus symbols: Error values e = \A n — A^ A \ for 
the sequence A n without Pade approximation. Solid lines and squares: Error values £ = \A^ A — Ai$~\ 
for the Pade approximant A„ A to the periodic orbit sums. 



d = 6, semiclassical resonances were calculated by the cycle expansion technique, including 
periodic orbits up to cycle length n = 13 || 0. In order to demonstrate the power of the 
PA we first apply it to the previously studied case with R = 1, d = 6. In billiards, which 
are scaling systems, the shape of periodic orbits does not depend on the energy E, and the 
classical action is given by the length L of the orbit (S po = hkL po ), where k = |k| = \l2MEjh 
is the absolute value of the wave vector to be quantized. 

We have calculated all periodic orbits with Maslov index /j po < 30, which corresponds 
to the set of orbits with cycle length n < 15. The resulting sequence of the partial sums 
{A n } of periodic orbit terms (eq. |]) converges for wave numbers k above the borderline 
Im k = —0.121557 Q which separates the domain of absolute convergence of the periodic 
orbit sum from the domain where analytic continuation is necessary, but strongly diverges 
deep in the complex plane, where the resonance poles are located. This is illustrated in fig. [l] 
for two different wave numbers k. The dashed line and the plus symbols in fig. [l]a show the 
convergence of the sequence {A n } at k = 150 — O.li. What is plotted is the error values 
e = \A n — Afjf'l, with the best known approximation to the true limit of the periodic orbit 
sum. As can be seen this sequence is slowly convergent, and about three significant digits are 
obtained at n = 15. The convergence can be accelerated using the PA, as is seen by the solid 
line and squares in fig.JJa showing the error values e — \A^ A — Af^\ for the sequence of the 
Pade approximants {A„ } to the periodic orbit sum. The periodic orbit sum has converged 
to six significant digits already by n — 9. The situation is much more dramatic in the deep 
complex plane (Im k < —0.122), e.g., at k = 150 — 0.5i. Here, the sequence of the partial sums 
{^4„} of periodic orbit terms exhibits exponential divergence, as can be seen by the dashed line 
and plus symbols in fig. [l]b. Nevertheless, this sequence converges when subjected to the PA 
implemented through the e-algorithm (see the solid line and squares in fig. |l|b). 

The resonances of the three-disk scattering systems have been obtained by a numerical 
two-dimensional root search in the complex fc-plane for the zeroes of the function l/g(k), 
where g(k) is the PA to the periodic orbit sum. A typical subset of the semiclassical resonances 
obtained is presented as crosses in fig. 0a. They agree perfectly with results obtained by the 



cycle expansion method in refs. p, M. For comparison the results of the cycle expansion |22 



are presented as squares. The dashed lines in fig. Ha denote the imaginary parts Im k = —0.1 
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Fig. 2. - (a) Semiclassical resonances (Ai subspace) for the three-disk scattering system with R = 1, 
d — 6. The crosses and squares mark the semiclassical resonances obtained by Pade approximant to 
the periodic orbit sum and by cycle expansion 1221 , respectively. The convergence properties of the 
periodic orbit sum and its Pade approximant arediscussed at Im k — —0.1 and Im k = —0.5 marked 
by the dashed lines (see text), (b) Same as (a) but with small distance d — 2.5 between the three 
disks. 



and Im k — —0.5, for which the convergence properties of the periodic orbit sum and the PA 
has been discussed above. 

Calculations similar to those discussed above for d = 6 have also been carried out at a 
shorter distance, d = 2.5, between the three disks. The semiclassical resonances obtained 
by PA to the periodic orbit sum and by cycle expansion |2^] are presented as crosses and 
squares respectively in fig. ||b. Again, the agreement is excellent with the exception of minor 
discrepancies for the imaginary parts of two resonances with Re k rj 94.3. It should be 
noted that the rate of convergence of the PA decreases when the distance d between the disks 
is reduced as is also the case when cycle-expansion methods are applied. A comparison of 
the convergence properties of the PA and the cycle-expansion for this specific system will be 
published elsewhere p3| . 

The circle billiard. - We now demonstrate the PA to periodic orbit sums for the example 
of the circle billiard. In sharp contrast to the three-disk system analyzed above, the circle 
billiard is an integrable and bound system, and, to the best of our knowledge, has not yet 
been treated by the cycle expansion technique Q or pseudo-orbit expansion The exact 
quantum mechanical eigenvalues E — Ti 2 k 2 /2M of the circle billiard are given by zeroes of 
Bessel functions J\ m \{kR) = 0, where m = 0, ±1, ±2, , ... is the angular momentum quantum 
number and R is the radius of the circle. The semiclassical eigenvalues can be obtained by an 
Einstein-Brillouin-Keller (EBK) torus quantization |24| resulting in the quantization condition 

m/kR) 2 — Iml arccos — - = n ( n H — ^ (5) 
kR \ 4/ 

where n = 0, 1, 2, ... is the radial quantum number. In the following we choose R = 1. The 
periodic orbits of the circle billiard are those orbits for which the angle between two bounces 
is a rational multiple of 2ir, i.e., the periods L po are obtained from the condition 

L po = 2m r sin7 , (6) 

with 7 = 7TTO0 jm r . Here, = 1, 2, ... is the number of turns of the orbit around the centre of 
the circle and m r = 2rn0, 2m ( j > + 1, ... is the number of reflections at the boundary of the circle. 
Periodic orbits with m r ^ 2m,/, can be traversed in two directions and thus have multiplicity 
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Fig. 3. - Real part (solid line) and imaginary part (dashed line) of the function l/g(k) for the circle 
billiard with radius R = 1 obtained by Pade approximant to the periodic orbit sum. The zeroes agree 
perfectly with the exact positions of the semiclassical eigenvalues (from eq. ph marked by the squares. 



2. For the amplitudes A po of the circle billiard, the Berry- Tabor formula || for integrable 
systems yields 

„4 po = ^i^/2L% 2 /m 2 r , (7) 

and the Maslov index is 

Upo = 3m r . (8) 

It follows from here that the partial periodic orbit sums A n in eq. (Q) should be calculated 
with order n = m r . (The choice n — fJ, po = 3m r is also possible but the numerical calculations 
would be less robust because the order of the polynomials in the PA are increased by a factor 
of 3.) We included all periodic orbits with m r < 100 in the PA to the function g(k). The real 
and imaginary parts of 1/ g(k) are presented as solid and dashed lines respectively in fig. |3|. The 
zeroes of the function \/g{k) agree perfectly to at least seven significant digits with the exact 
positions of the semiclassical eigenvalues obtained from eq. (]|) and marked by the squares in 
fig. | 

For chaotic systems the efficiency of semiclassical periodic orbit quantization can be es- 
timated by comparing the length of the longest orbit included in the calculation with the 
Heisenberg period Ljj = 2np given by the average density of states. However, for integrable 
systems the nearest neighbour spacing statistics is a Poisson distribution with high probability 
of nearly degenerate levels. The lengths of periodic orbits should therefore be compared with 
the length L max = 27r/Afc m i n (with Afc m ; n the smallest level spacing) rather than with the 
Heisenberg length. For the circle billiard, using orbits with period L po < 200, it is no 
problem to resolve level spacings Ak < 0.01, i.e., spacings where the characteristic length 
L max = 2n/Ak w 628 is by more than a factor of three larger than the maximum length of 
the periodic orbits. 

In conclusion, we have introduced the Pade approximant to perform summation of periodic 
orbit terms as a simple but powerful method for periodic orbit quantization. The Pade 
approximant allows the re-summation of the typically exponentially divergent terms of the 
semiclassical trace formulae. The method has been demonstrated on two systems with com- 
pletely different classical dynamics, viz. the classically chaotic three-disk scattering problem 
and the integrable circle billiard. The Pade approximant can be applied when the total periodic 
orbit sum can be divided into partial sums with respect to an integer ordering parameter n. 
It is an interesting speculation if the latter might be related to the cycle length of periodic 
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orbits if a symbolic dynamics exists, or, in the general case, to the Maslov indices of orbits. 
Clearly, an in-depth investigation of this point will be worthwhile but goes beyond the scope 
of this Letter. 

We thank G. Wunner and B. Eckhardt for fruitful discussions and comments on the 
manuscript. Financial support from the National Science Foundation (NSF-Grant No. PHY- 
9802534) is gratefully acknowledged. JM thanks HST for the kind hospitality during his visit 
at the University of Southern California where this work was initiated. 
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